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ABSTRACT 

We present an approximation for the average density profile of dark matter haloes 
in the ACDM cosmological model, which is accurate to within 10-15% even for large 
radial distances from 0.05i?vir up to lOiivir for halo masses ranging from 10^^'^ to 
j^Qis.o /j-ijv/g^ We propose a modified form of the Navarro, Frenk & White (NFW) 
approximation: p{r) = pnfwI?") + Air/Rvir)~^ + B{1 + r/i?vii)~^- This generalized 
expression, which is applicable to the external regions of dark matter haloes, only very 
slightly affects the density in the inner regions of haloes. The strong correlation among 
the different parameters in the model allows us to describe the profile in terms of just 
one parameter: the virial mass. We integrate our density profile to derive the enclosed 
mass in a sphere of a given radius and compare it with the NFW results. We find that 
the NFW underestimates the enclosed mass by more than 50% at 10i?vir, whereas our 
model reproduces the results from numerical simulations to within 2% accuracy even 
at this distance. We also use this new approximation to study the weak gravitational 
lensing and to obtain an analytic expression for the tangential shear. This allows us 
to quantify the contribution to the shear from the outer regions of the density profile. 
For the first time we calculate the difference between the tangential shear calculated 
via the NFW profile and the corresponding result when the external regions of haloes 
in cosmological simulations are taken into account. We find a 4% difference for all the 
mass ranges under study. 

Key words: dark matter - galaxies: haloes - large-scale structure of Universe - 
cosmology: theory - methods: iV-body simulations - methods: statistical 



1 INTRODUCTION 

Measurements of the mass distribution of dark matter 
(DM) haloes are essential for the theory of structure for- 
mation. This issue has been addressed in two ways: obser- 
vational and theoretical. The observational measurements 
have profited from recent gravitational lensing studies and 
also from measurements of galactic dy namics using large 
scale and massive galax y surveys (e.g. iMandelbaum et all 
120061 . IPrada et all '20031. Theoretical estimates made sub- 
stantial progress thanks to recent advances in both top- 
hat gravitational collapse models and fully non-linear evo- 
lution via A'^-body simulations (e.g. Sanchez-Conde et al.l 
12001 lAscasibar etaLll2007l , iDiemand et al.ll2007t ). 
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Cosmological A''-body simulations have been extensively 
used in order to provide predictions for the structure of DM 
haloes in the hierarchical clustering scenario. The Navarro, 
Frenk & White (NFW) density profile has become a con- 
venient formula for the description of DM haloes in a 
broad mass range. This profile provides a reasonable fit 
for the halo density distribution, despite that many other 
funct ions have be en suggested in order to impro ve accu- 
racy (iMoore et allflOTa. .Sersi3ll968l . lEinastd[l965l profiles. 



see also lMerritt et al.ll2005l for reference) . Nevertheless, the 
NFW density profile f N avarro et al.lll997h was proposed and 
contrasted to fit the inner regions of DM haloes roughly up 
to the vi rial radius, where h aloes are expected to be in equi- 
librium (ICole fc Lacevlll996l ). Therefore, it was not expected 
that the NFW analytical profile will describe the external 
density profile properly. Our aim here is to extend the va- 
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lidity of this approximation to describe t he outer regions of 
DM haloes. Moreover, recent works (see iPrada et alj |2006| . 
ICuesta et alj l2007l . IPiemand et alJ l2007h have shown that 
the DM virialized regions extend well beyond the formal 
virial radius especially for low-mass haloes. This reinforces 
our interest in studying haloes beyond the virial radius. 

On the other hand, current observational techniques al- 
low us to measure the mass distribution around galaxies 
and clusters at large distances. The dynamics of satellite 
galaxies (e.g. IZaritskv fc White! [l994l: IZaritskv et al.l[T997l : 



m 



IPrada et al.Tl2003l : lBrainerdll2004l : IConrov et al.ll200Vl ) pro- 
vide strong constraints on the shape of the density profiles 
of DM haloes in the outer regions. Another technique which 
has proven its usefulness despite its recent development, is 
the analysis of the weak lerising effect in gravitational lens- 
ing theory (e.g. Mellier 199ff; Bartelmann & Schneider 2 00 ll: 
Hockstra et al. 2004 Hcymans ct al. 2006; Masscy ct al.l 



20071 '). In both observational techniques, the external regions 



under study are often beyond the virial radius. Yet, the use 
of the NFW profile to model observations is still very com- 
mon in the determination of DM halo mass. Although the 
NFW profile is a great improvement as compared to the 
isothermal profile to describe for example th e weak lensing 
observations (e.g. iMandelbaum et al.ll2006a ). it is not suit- 
able for extrapolation at large distances beyond the virial ra- 
dius. The study of the physical processes occurring in the ex- 
ternal regions of DM haloes will provide unique insight into 
our understanding of the formation of DM h aloes. A good 
example is the recent observational work by iGeorge et al.l 
120081 who measured temperature and density profiles be- 
yond the virial radius of a cluster and this allowed them to 
produce improved constraints on the mass and gas fraction 
profiles. A first step in this direction is the study of aver- 
age density profiles in the outer regions of DM haloes in 
cosmological simulations. 

In this paper, we present an approximation for the den- 
sity profile of DM haloes well beyond the virial radius, ob- 
tained from high resolution cosmological simulations. This 
approximation focuses on the external regions of the haloes 
as an extension of the NFW formula that fits well the inner 
regions. We took the NFW model as our starting point be- 
cause of its simplicity and with the aim to derive an expres- 
sion which is easy to implement in many of the applications 
mentioned above. In order to describe representative DM 
halo density profiles and compare them with our approx- 
imation, we average over many individual density profiles 
for a given mass range. Hence, our conclusions apply for 
the ensemble of haloes of a given mass, with a halo-to-halo 
mass depende nt scatter. As it is a lready known from the 
early work of iNavarro et al.l l| 19961 ). individual haloes at a 
given range of mass usually show a similar profile. Thus, 
the analysis of average halo density profiles happens to be a 
reasonable option for our purposes. 

The approximation presented in this paper includes free 
parameters which have to be fixed for a given mass range. 
However, we will also study some correlations among them 
and prove that only one parameter (for example the virial 
mass) is enough in order to provide a high quality fitting 
function for the inner and outer regions. Moreover, our ap- 
proximation can be useful to obtain a suitable formula to 
describe the projected mass density of average DM haloes. 
We determine with this formula the tangential shear asso- 




Figure 1. Average matter density profile for a Milky Way size 
DM halo. The numerical profile given by cosmological simulations 
(dotted line) is compared to the fit to the NFW model (Eq. |2] 
long dashed line) , the mean matter density of the Universe (Eq. [S] 
short dashed line), and our model ('Eg. 1121 solid line). The upper 
panel shows the relative error for each approximation with respect 
to the numerical profile. 



dated to gravitational lensing. We also have compared our 
results with those obtained using the NFW formula, in order 
to assess its validity at large distances from the halo centre. 
The paper is organized as follows. In Section 2, we intro- 
duce our approximation that describes the density profiles of 
DM haloes well beyond i?vir. In Section 3, we fit the numer- 
ical density profiles drawn from cosmological simulations. 
In Section 4, we show the application of our model to the 
measurement of the tangential shear in gravitational lensing. 
Finally, we discuss and summarize the main conclusions in 
Section 5. 



2 THE THEORETICAL FRAMEWORK 

Let us start with a brief review of the NFW approximation 
for DM haloes in cosmolog ical A^-body simulations (see also 
e.g. iLokas fc Mamonll200ll ). This is a simple function which 
provides a good description for the density profile inside the 
virial radius. Yet, it turns out to give wrong results for the 
external regions since it does not take into account the mean 
matter density of the Universe, which is the m ain contribu- 
tion a t large distances from the halo centre (see lPrada et al.l 
I2OO6I ). The analytic expression for the NFW density profile 



PNFw(r) 



(1) 



{r/vs) (1 +r/rs) 

This formula describes the internal profile region quite 
well (up to one virial radius approximately) and depends on 
two parameters: the characteristic density p^^^ and the 
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characteristic scale radius 



The latter is the radius in 



which r^ -pNFwir) reaches its local maximum, i.e. ^"^^' ' = 
—2, and the former is just given by p^^^ — 4pNFw(rs). 
The NFW profile may also be described in terms of the 
virial mass Afvir = M{< -Rvir), and the concentration pa- 
rameter c = Rvir/rs- The quantity 7?vir is the virial radius 
of the halo, defined as the radius of a sphere enclosing a 
given overdensity, the value adopted here is A = 340. Hence, 
the relation between the virial radius and the virial mass is 
Mvir = ^3A0 p Rli^. The NFW profile can now be written 
as follows: 



pNPw(s) ~ 



(2) 



cs (1 + cs) 

where s is the radial coordinate scaled to -Rvir, i.e. s = 
r/RviT- Note that p^^^ is now a function (in general) of both 
c and Mvir. In order to find this function, let us calculate 
the halo enclosed mass from the density NFW profile: 



Mnfw(s) = 



3Mvir Ps 
340p t 



ln(cs + 1) - 



CS + 1 



(3) 



Since Mvir is defined as the mass inside the virial radius, 
we have. 



M{S ^1)= Mvir 



(4) 



It is straightforward to derive the expression for the 
characteristic density in terms of c and Mvir. It turns out 
that p^^^ is actually a function of the concentration pa- 
rameter only, i.e.: 



NFW 340 3 / \ - 
where the function g is defined as follows: 



ln(a; + 1) - x/{x + 1) 



(5) 



(6) 



Therefore, the mass inside a sphere of a given radius is com- 
pletely given by the virial mass of the halo and its concen- 
tration: 



M(s) 

Mvir 



5(c) 



ln(cs + 1) 



cs + 1 



9{c) 

g{cs) 



(7) 



To describe the external regions of DM haloes, it is es- 
sential to decide the kind of function that must be added to 
the NFW approximation, in such a way that its contribu- 
tion can be neglected in the in ner parts. A first a ttempt to 
improve this fit is proposed in IPrada et al.l l|2006l ). At very 
large radii the DM density profile should not tend to zero 
as the NFW profile does, but instead it should tend to the 
mean matter density of the Universe. Then, the modified 
profile is: 



Pmod(,S) = 



cs(i + csy 



+ P 



(8) 



The parameter p™° is no longer given by Eq. (5) . If we use 
the condition (4) we obtain in this case: 



mod iio3/\— ^ 

ps = 113 c g[c)p = ps 



nC 9{c)p 



(9) 



In addition, cosmological DM halo density profiles show 
a tiny or null influence on the concentration parameter at 
large distances (r > 7?vir), so that the density proflle in the 
halo outskirts can only depend on the remaining parameter 



Mvir. Therefore, in order to approximate the transition to 
the outer regions we add a function / which depends only 
on s (i.e. on the virial radius, but not on the concentra- 
tion). This makes Mvir the only relevant parameter at large 
distances from the halo centre. The new density profile ap- 
proximation that we propose here is a simple extension of 
the NFW formula, i.e., 

Ps 



P{s) 



cs (1 + cs)' 



+ P + f{s)p 



(10) 



where f(s) is the new function that we are proposing in 
order to improve the fit to the external parts of DM haloes. 
We choose the following rational function: 

fei 62 



f{s) 



(11) 



! S + 1 

There are several reasons for this choice. On the one hand, 
we aim for a simple formula which allows an easy calculation 
for the tangential shear. On the other hand, we want to 
minimize the influence of the new added parameters (&i and 
62) on the internal profile, where the NFW profile is steeper, 
except for the very inner regions in which NFW and f{s) are 
both proportional to r~^. We thus propose in this work the 
following model as an approximation of the density profile 
from the inner regions up to large distances well beyond 

Rvii'- 



p{s)^ 



+ 



7 + ;ti + '|^ 



(12) 



CS (1 + cs) 

It is worth noting that the enclosed mass and the character- 
istic density are now given by these relations: 



M{s) 



3 Mvir 

340 



^g{cs)p 



+ bi— + 



+ 62 ln(s -M) - s -f ■ 



1 1 


/ 


1 


77 + -bi -h 


In';- 




3 2 


^ 


2 



+ 



&2 



c g{c)p 



(13) 



(14) 



Applying the condition (J4]) to the expression (|13p leads 
to Eq. (I14|) . which quantifies the infiuence of external param- 
eters &i and 62 on the value of p^. In the case of &i = &2 = 0, 
which is equivalent to f{s) — 0, we recover the modified 
NFW characteristic density pT""^. 

An example of this new fitting formula for Milky Way 
size haloes is shown in Figure 1. Here we compare differ- 
ent approximations to flt the numerical density proflle. In 
order to reduce the range of variation along the y-axis, we 
plot the function (p/p) ■ (r/Rvir)^ instead of just p/p. The 
relative errors for each one of them are also displayed. The 
numerical profile given by cosmological simulations (dotted 
line) is compared to the NFW model (Eq. [2]), the modified 
NFW model (NFW plus the mean matter density of the Uni- 
verse, Eq. |5J and our model (Eq. I12|l . The fit in the inner 
part is quite similar for all the models as expected. There is 
an obvious discrepancy beyond the virial radius between the 
NFW flt, the modifled NFW approach, and the data from 
numerical simulations. Our model helps to alleviate this dis- 
agreement, providing a reasonable flt going from the inner 
regions to well beyond the virial radius up to lOflvir- 

The description of the DM distribution from small to 
large distances from the centre of a halo presents an added 
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Table 1. Results of the fitting of our approximation to average halo density profiles in different mass bins. The first column sets the name 
of the bin, the second shows the number of distinct haloes in each mass range, the third corresponds to the mass range in logarithmic 
scale, the fourth is the mean virial mass in this range, the fifth the concentration in the NFW model, the sixth, seventh and eighth show 
the best-fitting values of the parameters of our approximation. The ninth and tenth columns indicate the position of the local extrema 
in the numerical density profiles. 
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Figure 2. Halo-dark matter correlation function ^Halo-DM for 
average haloes in four different mass bins. The solid line represents 
our approximation whereas the dotted curves are the results from 
numerical simulations. 



relevance as it is related to the halo-da rk matter correla- 
tion function ^Haio-DM = '■P'-'^l'-p (e.g. iHavashi fc Whitel 
I2OO3). Therefore, our approximation for the density profiles 
in the outer regions of DM haloes is also showing informa- 
tion about this correlation function. In Figure 2 we display 
Cnaio-DM for different bins in the mass range covered in this 
paper (see Section 3) . In this plot we can clearly distinguish 
between the one-halo term, which corresponds to the halo 
density profile, and the two-halo term, which corresponds to 
the distribution of matter in its outskirts. Our approxima- 
tion provides a reasonable description in both regions. 



3 FIT TO SIMULATED DENSITY PROFILES 

To carry out our analysis we have used several high- 
resolution cosmological simulations. The initial conditions 



of these simulations have been set up from the power spec- 
trum corresponding to the ACDM model with the following 
values of the cosmological parameters: ^m— 0.3, ^a~ 0.7, h 
= 0.7 and as= 0.9. We used three different boxes whose size 
is 80, 120 and 250 h~^Mpc. The mass resolution in these 
boxes is 3.14 x 10*, 1.06 x 10^ and 9.59 x lO^fe'^M© respec- 
tively, and the space resolution is 1.2, 1.8 and 7.6 h~^'kpc, 
respectively. The final output of the simulation is obtained 
by tracking the evolution of 512'' p articles in each box t o 
z = 0, using the iV-body code ART (iKravtsov et al.lll997l '). 
The collapsed DM haloes are detected using Bound Density 
Maxima halofinder (BDM). 

We use mean halo density profiles averaged over seven 
different virial mass bins ranging from 5.35 x 10^^ h~^MQ 
to 4.47 X 10"ft"^MQ. There is a total of 17,220 haloes in 
all simulation boxes but we constrain this sample with the 
following selection criterion in order to select distinct haloes: 
the centre of all selected haloes must be further than the 
virial radius of every halo with higher mass. This criterion 
reduces the number of haloes to 16,679 (see ICuesta et alj 
I2OO7I for details on this halo sample). 

Figure 3 shows the best fit provided by our approxima- 
tion to some average density profiles from simulations. In 
this plot we can see that there are deviations of the order 
of 15-20% at the position Smax of the local maximum of 
p{r) ■ r^ . This deviation is inherited from the NFW model, 
where the fit to numerical density profiles usually show a 
noticeable disagreement at r = Ts. However, the transition 
from the internal to external regions at 2-37?vir of the halo 
displays a similar difference due to the fact that our simple 
model is not able to reproduce the steepest region in the 
density profile. In Table 1 we give the results of the best fit 
of our model to the simulated average halo density profile 
for the seven different mass bins. 

It is interesting to keep in mind the halo-to-halo varia- 
ti on of DM den s ity pr ofiles. This has been shown for example 
in IPrada et al.l l|2006l ) . In Figure 4 we show the root mean 
square in p{r) normalized to the density profile itself. Al- 
though only small variations (at the level of 20-30%) are 
present inside _Rvir among different haloes, the density pro- 
files show fluctuations larger than twice the density profile 
itself beyond 2_Rvir. Haloes corresponding to small masses 
present a larger scatter. 

A regression analysis between the parameters in our 
approximation allows us to find that there are strong cor- 
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Figure 3. Average density profiles in four different halo mass ranges. The solid line represents the best-fit of our approximation to the 
numerical density profile (dotted line). The deviation around r = rs is inherited to the NFW profile, which is about the same amount of 
the deviation at the transition to the outer regions of the halo at 2-3-Rvir- 



relations among them. These results are shown in Table 2. 
An interesting correlation occurs between the NFW con- 
centration (hereafter cnfw) and the parameter c in our 
model. This tight relation motivates the search for an in- 
terpretation of this parameter. We recall that in the NFW 
model the concentration is related to the radius in which 
the function p{r) ■ r presents its maximum, i.e. Smax = 
1/cnfw- It is tempting to explore whether the parameter 
c in our model admits a similar interpretation. The condi- 
tion d/ds (s^ ■ p/p) = is a sixth degree polynomial equa- 
tion with two real and positive solutions, Smax and Smin, 



corresponding to the local maximum and minimum of this 
function respectively. The numerical solutions for Smax and 
Smin are shown in the last two columns of Table 1. The cor- 
relation between cnfw = 1/smax and c turns out to be: 

1 



= 0.967(±0.004) ■ c + 0.53(±0.03) (15) 

.5 in ax 

with a Pearson correlation coefflcient of 1 — r^ = 7.3 x 10~^ 
showing this evident connection. This means that our pa- 
rameter c, which is no longer related to the position of the 
local maximum of p{r)r'^ due to the addition of /(s), still 



6 H. Tavio et al. 

Table 2. Results from the regression analysis. The first two columns display the name of the correlated variables (dependent and 
independent variable respectively). In the remaining columns, the value of Pearson's correlation coefficient, the coefficients of the regression 
equations and their standard deviations are shown. In the first four rows a linear regression is performed between both variables 



: mx + n), whereas in the last two rows a quadratic relationship is shown {y = px^ 



-n). 



y 



CNFW c 0.9985 0.66 0.15 

log(cNFw) log(Mvir) 0.9857 2.44 0.08 

log c log(Afvir) 0.9788 2.50 O.fl 

fei b2 0.9994 2.9 0.4 

62 c 0.9904 121 8 
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Figure 4. Halo-to-halo variation of DM density profiles in four 
different mass bins. Whereas the fluctuations compared to the 
magnitude of the density profile are small below l-Rvin they be- 
come very large around r > 2i?vir- At large distances, this scatter 
is higher for decreasing halo mass. This effect is not an artifact of 
the different behaviour of mean density profile at large distances, 
but the intrinsic scatter due to halo environment. 



retains its original interpretation, with only a small varia- 
tion in its value. The function f{s) has therefore very little 
influence in the inner regions, i.e. in the NFW term of our 
approximation. 

The correlation between the NFW conc entration and 
the v irial mass is already well-known (e.g. iMaccio et al.l 
[20071). Therefore, the relationship between our parameter 
c and Mvir is straightforward as shown in Figure 5 for our 
sample average halo profiles. This correlation is extremely 
useful as it will allow us to derive simple mathematical ex- 
pressions for the outer density profile in terms of the param- 
eter c, which would be more complicated if they are written 
in terms of the virial mass. However, it is important to re- 
mark that this choice is just for our convenience, as we know 
that the parameter c has no effect in the outer halo regions. 



Figure 5. The correlation between the average virial mass in 
the mass bins under study and the concentration parameters, i.e. 
both the usual NFW concentration (solid circles) and the c in our 
approximation (open circles). Solid and dashed lines represent the 
best fit to a power-law for CNPw(-'^vir a^nd c(Mvir, respectively. 



Something more surprising arises when we examine the 
relation among the fitted parameters 61, 62 and c (as it can 
be readily seen from the last two rows in Table 2). This 
allows us to remove two parameters from our density model 
(Eq. 12) as &i and 62 can be expressed in terms of c with 
little scatter. 

Thus, taking into account that the parameters ps, bi 
and 62 are a function of c, expression (I12|l turns into: 

P(£) 
P 



ps{c)/p ^ bi{c) ^ b2{c) _^ ^ 



cs {cs + 1)2 ' s ' s + 1 
ps/p (fei + &2)s + bi 



cs (cs + 1)2 



+ ■ 



s(s+l) 



+ 1 



(16) 



where ps{c)/p, &i(c) and &2(c) are no longer free parameters 
but fixed in terms of c. This dependence is determined by 
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^ 40 



Figure 6. Top panel: correlation between the parameters fei and 
62 in our approximation. There is a tight linear correlation be- 
tween both parameters which describe the outskirts of the density 
profile, so our model gains simplicity. Bottom panel: correlation 
between the parameters 62 and c in our approximation. Although 
a linear regression fails in order to find a correlation of this two 
parameters, the quadratic regression works reasonably well. 



the expression (|14p and the values shown in the last two 
rows in Table 2, i.e., 



Ps/p{c) = 
62(C) = 



61(c) 



0.167c^ 
0.579 c^ 



= -0.558 c^ 



-4.12 c +143.4 
-14.94 c +120.8 
+14.02 c -107.4 



Table 3. Results of the constrained fitting of our approximation 
to halo density profiles. As opposed to Table 1, here c is the only 
free parameter, whereas the values of 61(c) and 62 (c) are obtained 
using Eq. I I17I I for the best-fit value of c. We include the fourth 
and fifth column to allow direct comparison with Table 1. 



Name 


(Me/h) 


c 


fei(c) 


b2{c) 


HI 


5.35 X 10" 


12.42 


-19.32 


24.60 


H2 


1.64 X 10^2 


10.12 


-22.66 


28.95 


H3 


5.08 X 10^2 


8.77 


-27.37 


34.35 


H4 


1.58 X I013 


7.12 


-35.89 


43.82 


H5 


4.82 X 10^3 


6.80 


-37.89 


46.02 


H6 


1.48 X 1014 


6.17 


-42.17 


50.70 


H7 


4.47 X 10" 


5.26 


-49.12 


58.27 



It is important to realize that Eq. (|16|l has been grouped in 
such a way that each term is related to a region of the halo 
density profile (internal, middle and external, respectively). 
After the elimination of 61 and 62, we obtain again the 
best-fit for the only free parameter c, taking into account 
that now the values of &i and 62 are fixed and can be ob- 
tained using Eq. (|17p . We note that Eq. (|17|l implicitly in- 
cludes the correlation between c and Mvir, as the use of 
the parameter c instead of Mvir makes this expression much 
simpler. We have chosen this form for our mere convenience 
although c has nothing to do with the outer halo regions. 
The results are given in Table 3, where the fourth and fifth 
column allows a comparison with Table 1. Only small dif- 
ferences about few per cent are present between the best-fit 
values for 61 and 62 in Table 1 and the new results for the 
best-fit value of c in Table 3. 

Apart from the discrepancies already shown in Figure 3, 
we also note that the slope in the density profiles from the 
simulations tends to be shallower as compared to our ap- 
proximation starting around 8 virial radii. We therefore ad- 
dress the issue of the range of validity of our approximation. 
In order to study this, we will use a density profile which 
extends up to nearly 30 virial radii (see Fig. 7). 

At these very large distances we are able to analyse the 
asymptotic behaviour of both numerical and approximated 
density profiles. To this aim, it is more useful to use the func- 
tion y — (p/p) s since we can study its oblique asymptote. 
If we multiply Eq. (|12|) by s we obtain that the asymptotic 
behaviour is: 



y = s + bi + b2 



(18) 



(17) 



The leading term at large distances corresponds to the 
trend of the density profile to have the same value as the 
mean matter density in the Universe, but there is also a 
dependence on the coefficients 61 and 62. In Figure 8 we 
can also see that this dependence makes our approximation 
slightly shifted by a constant value 61 + 62 with respect to 
the data from cosmological simulations. This means that 
the 0{l/r) terms in our approximation are still important 
in the range from 10 to ~ 30-Rvir just before entering the 
asymptotic regime where the density profile tends to the 
constant mean matter density p. 

This shows the limitations of our approximation: while 
it provides a reasonable fit up to 10 virial radii, the goodness 
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Figure 7. The density profile of an average lialo witfi Mvir — 
5x 10^^ h~^MQ up to nearly 30i?vir- Although our approximation 
implements the right trend at large distances from the halo centre, 
residuals of the order of 1/r are amplified since the density profile 
is multiplied here by r^. 




Figure 8. The density profile p/p times r/iJvir versus r/iJvir 
for our approximation (solid line) and the data from numerical 
simulations (solid circles). The dashed and the dotted lines rep- 
resent the oblique asymptotes which correspond to our fit and 
the simulation data, respectively. The vertical shift between both 
asymptotes shows the presence of 1/r errors in our approxima- 
tion. 



Figure 9. Enclosed mass in a sphere of radius r versus r/iJvir 
for the average halo with Mvir — W^^h~^MQ. The dotted curve 
represents the simulation data, long dashed and short dashed lines 
show the plain NFW and the modified NFW profile (i.e. NFW 
plus the mean matter density p) respectively. The solid line is 
our approximation. Top panel represents the relative errors with 
respect to the simulation data. 



of the fit is not so good at larger distances. Tliis is mainly 
due to the simplicity of the approximation, wliicli is one 
of the premises in our model. Of course, we can always add 
more rational terms to our function / but this would compli- 
cate the model and the interpretation of the different terms. 
However, the density at distances larger than lOiJvir is sim- 
ilar to the mean matter density, e.g. p{s > 8) < 1.5p. With 
this in mind, we can calculate the enclosed mass in a sphere 
of a given radius, and compare with our approximation (see 
Figure 9). We can see that even comparing with the mod- 
ified NFW density profile (i.e. NFW plus the mean matter 
density in the Universe), our approximation is a significant 
improvement regarding the enclosed mass. While the modi- 
fied NFW profile can overestimate the enclosed mass by 5% 
at 3-Rvir, and underestimate it by 10% at SRvh, our approx- 
imation keeps this uncertainty below 2% even at 10i?vir, 
which makes our simple model more than enough for this 
application. On the contrary, the plain NFW profile is not 
suitable beyond 4_Rvir in order to estimate the enclosed mass, 
as the estimated mass falls short by more than a factor of 
two with respect to the data from the simulations at 10J?vir 
(see Figure 9). 



4 APPLICATION TO GRAVITATIONAL 

LENSING 

One of the most relevant experimental techniques for deter- 
mining the distribution of dark matter is based on the gravi- 
tational lensing effect (for some specific applications, see e.g. 
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Gavazzi et al.1 20031. Broadhurst et al.1 20051. Johnston et alj 



20071 . iMandelJjaumet al.ll2006l . lMandelbaum et al.ll2008l ). In 

these observations the quantity in which we are interested 
is the tangential shear 74 , which describes the image elonga- 
tion perpendicular to the line connecting the image from the 
source and the distorting lens. This quantity is related with 
the surface mass density E (projected on the lens plane, 
which is perpendicular to the trajectory of the incoming 
ray), through the following formula: 



lt{R) 



AE _ E(< R) - T.{R) 



(19) 



where R is the radial coordinate projected on the plane of 
the lens (impact parameter). E(< R) is the mean surface 
density enclosed in a radius R, and E^rit is the critical sur- 
face density, which depends from the distance between the 
lens and the source, i.e., 



Ds 



47tG DlDls 



(20) 



where Ds, Dl and Dls are the distance to the source, the 
distance to the lens, and the distance between the source 
and the lens, respectively. 

In fact, the quantity which can be observed is the shear 
in units of critical density AE(_R). We can derive its ex- 
pression using our approximation for the dark matter halo 
density profile. 

The surface mass density is given by the following for- 
mula, i.e. the projection of the volume density p{r) in the 
line of sight: 



E(i?) = 2 /' 

J R 



p{r) 



\fr 



7?2 



dr. 



(21) 



instead, the mean surface density is defined as follows: 

c-R 



E(< R) = ^ [ E(i?') dR'. 
R Jo 



(22) 



For simplicity in our calculations, we use distances which are 



scaled to the virial radius (s = -^^— ; ^ 
the last two equations transform into: 

n+00 



E(5') = 2R^ 



Pis) 



Vs 



52 



ds 



Therefore, 



(23) 



n<S) = ^j E{S')dS' 



(24) 



where we choose for p{s) our approximation for the den- 
sity profile given by Eq. (|16|) . in which 61 and 62 are func- 
tions which depend on our parameter c as explained in the 
previous section, i.e.. 



pis)^ 




&l(c) &2(C) 



(25) 



d(s) 



Now we compute the tangential shear for this density 
profile. As we can separate the contributions from differ- 
ent terms in Eq. (|25|l due to linearity of the integrals in 
the definition of the shear, we will calculate each term (i.e. 
the internal, the middle region and the contribution from 
background density) separately. The last term has associ- 
ated a trivial contribution, because it is a known fact that 
for any homogeneous density field the tangential shear is 



null. Therefore, in the case of the background density we 
obtain AE(,(5') = 0. 



4.1 Internal tangential shear 

The first term in this equation is formally identical to the 
NFW model and the surface ma ss density E for th is model 
has already been determined bv iBartelmannI (|l996r ): 



Eint('j) 



2psRvir 



ds 



(cs-(-l)Vs2 -52 



Eint(o) — < 

where (p is given by: 



c c^S-^ — 1 ' c 



(26) 

(27) 



3c 



ff s= ^ 



atanhvi^r^ „ ^ ^ , 
, ^ — < a: < 1 



{x)^l 



arctanVi2_i 



a; = 1 



X > 1 



(28) 



On the other hand, we get from the calculation of the mean 
surface density: 



Ei,it(< S) — 



2psRvir 



' Jo '-^P^dS' 0<S<l/c 

1/c S = 1/c (29) 



and hence we can write this in terms of 0, i.e.: 
E„,(< S) = ^P^-^ ^{cS). 



(30) 



Therefore the tangential shear AEi„t(S') — Ei„t(< S) 
Ei,it(o) is: 



AE„,t(5') = 



2ps_Rvir 






-^'^(-^^-^ S/l/c 



(31) 



2/3 S = 1/c 



So, if we introduce for our convenience the function 7/) 
defined as follows: 



■,p{x) 



X (h(x) — l I 1 



2/3 X = 1 



(32) 



We can simply write the tangential shear for the internal 
region in this form: 



AE„,,(5) = ^^^-5^ i,{cS) 



4.2 Middle tangential shear 



(33) 



The second and third terms in the equation (|25|l rep- 
resent the contribution from the intermediate region of 
the density profile to the tangential shear. For pmid{s) = 
( ^^ -1- -j^ ) p the surface density is: 
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^mid\^) — ^-rtvirP 



— 2i?virP 



(&i(c) + b2(e)) s ds 

S (s + l)Vs2-5'2 



(fei + ba) lim In ^^' /' + 1 + 0(g) 

5— . + CX) O 



where 0(3) is defined in the same way as above, and for 
the mean surface density we get: 



Smw(< S) — 2_RvirP 



+ (bi + &2) hm In 

^— +CXD 



(6l+&2) + 



5 



+ 11 0(5') dS' 



The tangential shear from the intermediate region is there- 
fore: 



^^mid{S) — 2RvirP 



&i+fe2(l+<?i(S)-| I <f){S')dS' 



(34) 



Unfortunately there is no primitive function for which 
could be written in terms of elementary functions. Never- 
theless, it is possible to find a fitting function for it. We will 
fit i /jf (j}{S')dS' - (j){S) over the interval s G [0, 10] using 



the function 



l + (a:/a) 



The result is: 



(t){S')dS' - (t>{S) 



1 



{£)' 



(35) 



We are now in a position to write down the complete 
expression of the tangential shear in units of the critical 
density for our model, i.e.. 



AS(S) = 2i?virP 



f-c^^'^'^^'^^'^ijirs 



(36) 



We can now compare this tangential shear with that 
one given by the NFW model (which is formally identical to 
AEi„t(5'), except that we have to replace ps and c by p^^^ 
and cnpw respectively): 



AE^ 



r{S) 



2pl 



Rvir 



cnfw 



■i/)(cNFwS') 



(37) 



In Figure 10 we show the relative difference between 
AEnewC^) (which does not account for external region of 
the halo) and AE(S') as given by our model for the differ- 
ent halo mass bins. The differences are small and always 
below 4%. They are not appreciable at the level of current 
experimental sensitivity. As expected, the contribution from 
the matter distribution in the outer regions of DM haloes is 
small, but for the first time it has been estimated to what 
extent it may have infiuence in the measurement of the tan- 
gential shear. Besides, we find an obvious trend with halo 
mass: the difference between both AEnfw(S) and AE(S) 
is higher for most massive haloes. 



5 DISCUSSION AND CONCLUSIONS 

In this paper we have presented a simple approximation 
for the DM density profiles of haloes with masses ranging 
from 10^^'' to W^'^-°h~'^MQ, which is valid even beyond the 
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Figure 10. The relative difference between the tangential shear 
(in units of critical surface density) for NFW approximation and 
our model, (ASnfw(S) - AE(5))/AS(5), in four different mass 
ranges. 



10i?vir 



s{l + cs) 



The expression we are sug- 
^ + ^ + l") -O' where 



+ 



virial radius up to 
gesting here is p{s) 

s = r/Rvir- This approximation is an extension of the NFW 
formula but includes two additional parameters. We have 
shown that these parameters are very well correlated with 
the virial mass, so that the density profile is effectively just a 
function of Mvir- Other approximations found in the litera- 
ture will prove to be useful to describe, with high accuracy, 
the halo density profile up to the virial radius. Yet, most 
of them fail when they are applied to fit the outer regions, 
where the mass predicted by extrapolation of these profiles 
is far below the actual mass in these regions. 

The description of DM distribution far away from the 
halo centre is especially interesting. In particular, the halo- 
dark matter correlation function is related in a straightfor- 
ward way to the average density profile. Although the tran- 
sition from the one-halo term to the two-halo term in this 
correlation function is present before 10i?vir, our approxi- 
mation has turned out to be a reasonable description of the 
DM distribution even at these distances. In order to build 
the numerical mean density profiles, we averaged over many 
hundreds of haloes from high-resolution cosmological sim- 
ulations so that the profiles corresponding to most of our 
mass bins are entirely unaffected by statistics. This proce- 
dure for averaging density profiles is analogue to the stack- 
ing meth od, used in obse rvational studies like the one by 
iMandelb aum et al.l (|2008l ) to infer the density profile of a 
cluster of galaxies. This similarity is useful for the compari- 
son of the results from cosmological simulations to the real 
data. 

This parametrization for the average density profile is 
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accurate to within 10-15% in the range from 0.05i?vir to 
107?vir- There are two main discrepancies from the numer- 
ical density profile which have a different origin: whereas 
the overestimation around r = r^ is inherited from NFW 
profile, the overestimation just beyond the virial radius sug- 
gests that our model is not able to reproduce the steep- 
est region. This steep region just outside J?vir is more pro- 
nounced for most massive haloes, suggesting a depletion of 
the halo outskirts due t o dark matter infall (|Prada et alj 
I2OO6I . ICuesta et al.ll2007l ). The presence of our additional 
terms with respect to the NFW formula has only a very 
small influence on the inner regions of the density profile, so 
that our approximation can also be considered as an exten- 
sion of the NFW profile. At larger distances our model shows 
deviations around 20% in the range 10-30-Rvir just before en- 
tering the asymptotic regime. These deviations are caused 
by our additional (r/_Rvir)~^ terms which improve the fit 
in the interesting region below 10-Rvir where the density is 
much higher. In any case, we must remark that our approx- 
imation implements the correct asymptotic behaviour: the 
density profile tends to the asymptotic value of the mean 
matter density of the Universe p. 

The cumulative mass inside a sphere of a given radius 
is underestimated by more than 50% at 10i?vir by the NFW 
formula. On the contrary, it is much better approximated 
(to within 12%) when the NFW profile is modified by addi- 
tion of the mean matter density, although with our model 
the difference with numerical density profiles is reduced even 
up to 1% in the range l-9i?vir. This is especially interesting 
for new measurements of the enclosed mass beyond virial 
radius in X-ray clusters IjGeorge et al.l I2OO8I ) , where plain 
NFW is still used even at r > 17?vir. While current obser- 
vations cannot distinguish between modified NFW and our 
approximation, in the near future they should be able to 
find the need for adding the mean matter density term to 
the density profile. 

We have also presented an application for our ap- 
proximation in the context of mass estimation using 
gravitational lensing effect. We derived expressions for 
tangential shear corresponding to different regions around 
the halo, which are in turn related to the different terms 
in our approximation. The contribution from the outer 
regions is small as compared to the contribution of the 
inner region, as expected. We calculated the difference 
between this tangential shear and the one derived from the 
NFW profile as a function of distance, showing that the 
inclusion of the outer regions produces a difference around 
4% . This small difference could provide an observational 
test for the validity of our approximation, which has 
been derived from the results of cosmological A'^-body 
simulations. Although present resolution of weak lensing 
experiments prevents us from drawing a robust conclu- 
sion, the stacking of different observations should prove 
that this approximation, which includes the contribution 
of external regions, is more realistic than most of the so 
far proposed density profiles, which do not account for them. 
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